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1 A path to loop quantum gravity 

In order to motivate the loop approach, it is useful to begin with a review of the two main 
paths to the quantization of general relativity : the covariant, or functional integral, approach, 
and the canonical, or Hamiltonian, approach. Their partial successes and difficulties will shed 
light on the problem of quantum gravity. 

1.1 Covariant approach 

The basic goal of the covariant approach is to define a functional integral for General 
Relativity, 

'"p5^,e-^^--(s), (1.1) 



where 

Seh{9^u) ^^^JA^xV~99''''Rt.A^{9)) (1-2) 

is the Einstein-Hilbert action. Here g is the determinant of the metric, and T{g) is the Levi- 
Civita connection entering the covariant derivative of vectors and tensors, 

V^v^^d^v^ + T;^{g)vP, r;;^(.9) = ^.9'^^[9pgAp+5,,gAp-aAgp^]. (1.3) 

From courses in QFT, we know how to properly define only the Gaussian integral 

Vipc^pf- f^ipDipj. (1.4) 

This corresponds to a free theory. An interacting theory can be treated in perturbation theory 
from a generating functional based on (1.4). However, there is no quadratic term in the 
Einstein-Hilbert action (in fact, it is not even polynomial). This obstruction can be bypas- 
sed if one performs a perturbative expansion of the metric. Consider the field redefinition 

9fj.u = Vt^iy + hfj.u, (1-5) 

where T/^y is fixed to be the Minkowski metric, and /i^^ the new dynamical field. The intro- 
duction by hand of the background field jy^^ is crucial : if we treat h^iy as a small fluctuation, 
a perturbative expansion of the action gives (in the De Donder gauge) 

Seh{9^.u) = Tp^ I d"^ K,nh^^'' + o(/i3). (1.6) 

The background metric can be also used to define the Wick rotation. This action is amenable 
to the machinery of QFT. Among the main results of this background-dependent perturbative 
approach, let us highlight the following ones : 

1 . The quanta of the gravitational field are interpreted as massless spin-2 particles (the 
"gravitons"). 
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In the low energy static limit, one can compute quantum corrections to the classical 
Newton potential [1], recovering the classical relativistic correction o{G^) and a new, 
purely quantum one o{hG^) : 

r?ll?7l2 / „'7T,l+?Tt2 127 hG _o, 

307r2 -2 



Vn = -G^-^[1 + G^ ±-—-— + oir-')] (1.7) 



3. As for other quantum field theories, one can define the S-matrix to describe the funda- 
mental observables, and compute scattering amplitudes of gravitons with matter fields 
and with themselves [2] 

This brief list does not exhaust the interesting results obtained with this approach. On 
the other hand, there is a fundamental difficulty with it. As for any quantum field theory, 
infinities appear when considering the effects of arbitrarily small ("ultraviolet") field fluc- 
tuations. These divergences are usually dealt with via the procedure of renormalization. This 
technique, successful for X(j)'^ or gauge theories, fails in the case of gravity. This can be ex- 
pected from dimensional arguments, and the rigorous proof that the perturbative quantization 
of general relativity fails because of non-renormalizable ultraviolet divergences was obtained 
in the late eighties, by Goroff and Sagnotti [3]. 

This means that the theory constructed as sketched above can be used for low-energy 
calculations, but it would be inconsistent if taken seriously at all energy scales. Heuristically, 
the inconsistence can be explained as follows. The key field redefinition (1 .5) assumes that we 
can quantize the field /i^^, and that its dynamics takes place on a fixed classical background 
Ty^i/. However, as we increase the energy of the fluctuations its backreaction increases, and 
it becomes inconsistent to assume that the background stays fixed, unperturbed. Although 
gravitons capture correctly the low-energy physics of the gravitational field, they might not 
be the right quantities to describe the quanta of gravity at Planckian energies. 



1.2 Canonical approach 

The canonical formalism is based on second quantization, and describes the quantum 
theory in terms of functionals of the fields, e.g. '!'[(/)] for a scalar field 0, like the familiar 
function of the configuration variables ip[x] in quantum mechanics. The dynamics is descri- 
bed by the quantum Hamiltonian H and the Schodinger equation 

*;i^vi/[0] = ij L A") vi/[0] (1.8) 

For the reader less familiar with this approach, see [4]. Given the overwhelming success 
of the functional integral approach to the quantization of (non-gravitational) field theories, 
the canonical formalism is often underappreciated in many courses, and a certain taste of 
"old fashion" is associated with it. It was applied to general relativity by Arnowitt, Deser 
and Misner (ADM), Dirac, Wheeler and De Witt, among many others. The loop approach is 
related to this original idea, which we will now review in some details. For more on it, see 
e.g. [5]. 
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1.2.1 ADM formalism 

In order to put the Einstein-Hilbert action into canonical form one needs to identify the 
variables which are canonically conjugated, and then perform the Legendre transform. To 
that end, we make the assumption that M has the topology A^ ^ M x S where S is a 
fixed three-dimensional manifold of arbitrary topology and spacelike signature. This poses 
no restrictions if we require that M has no causally disconnected region (more precisely that 
M is globally hyperbolic). In fact by a theorem due to Geroch and improved by Bernal and 
Sanchez [6, 7] if the space-time is globally hyperbolic then it is necessarily of this kind of 
topology. 

Having made this assumption, one knows that Jv[ foliates into a one-parameter family 
of hypersurfaces St = ^ti"^) embeddings of E in A^. The foliation allows us to identify 
the coordinate t G M as a time parameter Notice however that this "time" should not be 
regarded as an absolute quantity, because of the diffeomorphism invariance of the action. A 
diffeomorphism g DifF [M.) maps a foliation X into a new one X' = X o (p, with a new 
time parameter t'. Conversely, we can write a general diffeomorphism (f> E Diff (Ai) as com- 
position of different foliations, (j) — X' o X^^. Hence, we can always work with a chosen 
foliation, but the diffeomorphism invariance of the theory will guarantee that physical quan- 
tities are independent of this choice. A theory of spacetime which unlike general relativity 
has a preferred foliation and thus a preferred time, is a theory which breaks diffeomorphism 
invariance. 

Given a foliation Xt and adapted ADM coordinates {t, x), we can define the time flow 
vector 

_ dXj^ix) 



'(x) 



dt 



(1,0,0,0). 



(1.9) 



This vector should not be confused with the unit normal vector to E, which we denote ri^. 
They are both timelike, q^uT^t'^ = 5oo ^nd g^v'n^ri'^ = ^1. but they are not parallel in 
general. Let us decompose r^ into its normal and tangential parts. 



'(x) =N{x)n^'{x)+N^'{x). 



(1.10) 




It is convenient to pai-ametrize n'^ = (1/iV, -N°-/N), so that N^^ = (0, N"-). N is called 
lapse function, and iV" shift vector In terms of lapse and shift, we have 



ff^.r'^T^ = .goo = -N' + gabN'^N", 



gOa = gabN'' = Na 
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Using these results, the metric tensor can be written as 

ds^ = g^^dx^'dx" = -{N^ - NaN") dt^ + 2Nadtdx'' + .gabda;"dx^ (1.11) 

where a ~ 1, 2, 3 are spatial indices and are contracted with the 3-dimensional metric gab- 

Notice that the spatial part gat is not in general the intrinsic metric on St. The latter is 
given by 

Qf^i^ ^ gf^iyn^n^. (1.12) 

Tensors on the spatial slide S, since their scalar product with n vanishes, can be equivalently 
contracted with g or q. The quantity g(^ = g^'^qpi, acts as a projector on St, allowing us to 
define the tensorial calculus on St from the one on M.. An important quantity is the extrinsic 
curvature of St, 

^M- = <'?£:' V^'n,'. (1.13) 

This tensor is symmetric and it is connected to the Lie derivative of the intrinsic metric 
^nQfii' = '^Kfj^iy. It enters the relation between the Riemann tensor of St (7?,) and that of 
M(R), 

Kp. = <l'X<l'p€'K''p'a' - K.aK>; - K^pK^ (1.14) 

This formula, proved in the Appendix and known as the Gauss-Codazzi equation, lets us 
rewrite the action ( 1 .2) in the following form, ' 

S = I dt I d^Xy/^N [7^ - K"^ + Tr{KK)] . (1.15) 

A key consequence of this analysis is the fact that time derivatives of N and N°- do not 
enter in the Lagrangian L. This implies that N and A^" are Lagrange multipliers, with null 
conjugate momenta SL/SN ~ 6L/5N°' = 0. The true dynamical variables are the spatial 
components qab only, with conjugate momenta 



^'^^ = ^ = Vg (A'"'' - Kq"'') . (1.16) 

oqab 
One can finally evaluate the Legendre transform, which after some calculations gives 

Seh {qab,7T''\N,N'') ^ ^J-^ f dt f d^'x [TT'^'qab-N'^Ha-NH], (1.17) 

where (the covariant derivative V contains only spatial indexes) 

Ha = -2^\/b(^) (1.18) 



and 



H — —p GabcdT^"" T^'^ - \^Tl, Gabcd = qacqbd + qadqbc - qabqcd- (1-19) 

V9 



1. We are assuming for simplicity of exposition that S lias no boundary, the case with boundary is treated 
extensively in [8] and [9]. 
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The quantity Gated is often called the supermetric, or DeWitt metric. The variation of the 
action (1.17) with respect to the Lagrange multipliers gives the equations 

Ha{q,7:)^0, H{q,7T)^0, (1.20) 

which are called respectively vector, or space-diffeomorphism constraint, and scalar or Ha- 
miltonian constraint. In the following we will also use the compact notation H^ = (H, Ha). 
Physical configurations, also called on-shell configurations with particle physics jargon, must 
satisfy these constraints. 

From ( 1 . 17) we see that the Hamiltonian of general relativity is 

Ii=T:^ fd'^xN'^Ha + NH. (1.21) 

IdttG J 

This Hamiltonian is peculiar, since it is proportional to the Lagrange multipliers and thus 
vanishes on-shell. Hence, there is no dynamics and no physical evolution in the time t. This 
puzzling absence of a physical Hamiltonian is in fact a consequence of what we discussed 
earlier ; the diffeomorphism-invariance of the theory tells us that t is a mere parameter devoid 
of an absolute physical meaning, thus there is no physical dynamics in t. This is the root of 
the problem of time in general relativity. For discussions of this problem, see [10] and [11]. 

1.2.2 Symplectic structure 

From now on, we work in units I&ttG — 1. The Hamiltonian formulation (1.17) allows 
us to study the phase space of general relativity. It is parametrized by the pair {qab,T^°'^), with 
canonical Poisson brackets 

{^''\t,x),q,d{t,x')] = SlS^^-^^ix - x'). (1.22) 

From there we can evaluate the following brackets among the constraints (for the explicit 
calculation see in the Appendix) 

{Ha{x),Ht{v)} = Ha{y)db5{x - y) - Ht{x)d'J{x - y) 
{Haix),H{y)} = H{x)daS{x - y) (1.23) 

{H{x),H{y)} = H'^{y)dJ{x ^ y) - H-{y)d'J{x ~ y) 

Notice that the right-hand sides vanish on the constraint surface (1.20). This means that the 
Poisson flows generated by the constraints preserve the constraint hypersurface. Constraints 
with this characteristic are said to ht first class, as opposed to second class constraints whose 
Poisson brackets do not vanish on-shell. First class constraints generate gauge transforma- 
tions on the constraint surface (See e.g. [12] for details). 

To see what the gauge transformations look like in our case, consider the smearing of the 
constraints 

H{N) = / H''{x)Na{x)d^x, H{N) = / H{x)N{x)d^x. (1.24) 
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An explicit computation (see Appendix) shows that 

[H{N),qab} - Cj^qab, {J?(7V), tt""} = Cj^tt'^' (1.25) 

which means that the vector constraint is the generator of space-diffeomorphism on S. The 
situation is somewhat subtler for the Hamiltonian constraint. We now have (see Appendix) 

{H{N),qab} = CnNqab, (1.26) 

{H{N),n-'} = A-jattt'"' + ig'-'^iViJ - 27VV99''"9''''^cd (1.27) 

The first bracket is the action of time diffeomorphisms on qab- The second bracket gives the 
action of time diffeomorphisms on Tiab, but contains also two extra pieces. These vanish if 
H = and Red = 0, namely on the constraint surface and for physical solutions (recall that 
in vacuum Einstein's equations read i?^^ = 0). Therefore, we conclude that the constraints 
H^ are the generators of the spacetime diffeomorphism group DifF(A^ ) on physical configu- 
rations. 

For general configurations, ( 1 .23) defines the algebra of hypersurface deformations, often 
called Dirac algebra or Bargmann-Komar algebra. A characteristic of this algebra is that it is 
not a Lie algebra. In fact, let us look at the smeared Poisson bracket 

{H{Ni),H{N2)} = H {g''\NidbN2 - N^dbN,)) . (1.28) 

We see that outside the constraint surface, the "structure constants" on the right-hand side 
contain the field g"-^ itself. Hence they are not constants at all, and the algebra (1.23) is 
not a Lie algebra, unlike DifF(A^). Instead, (1.23) shows that when we introduce a foliation 
to define the canonical formalism, we still have the symmetry of diffeomorphisms, which 
acts changing the foliation, and this new one, which acts deforming it. The two symmetries 
coincide on physical configurations. 

1.2.3 Constraints and physical degrees of freedom 

A simple counting of the number of degrees of freedom of general relativity can be done 
in the covariant perturbative approach described in Section 1 . L The structure of the linearized 
field equations shows that only two components of /i^^ propagate, which correspond to the 
two heUcities of a massless spin 2 particle. 

An advantage of the canonical formalism is that it allows us to confirm this counting in a 
more general and robust way. Recall in fact that in classical physics each point in phase space 
(i.e. "initial position and momentum") characterizes a physical trajectory, and the number of 
degrees of freedom is defined to be half the dimensionality of the phase space. In constrained 
theories, such as general relativity but also gauge theories, one has to be careful with the 
constraints. To that end, it is customary to distinguish a notion of kinematical phase space, 
and physical phase space. The kinematical phase space is the one defined by the Poisson 
structure of the theory. In our case, the space {qabi tt"^), with Poisson brackets (1.22). The 
dimension of this space is (6 + 6) • oo'^ = 12 . od-^. 

On this space, the constraints define a hypersurface where they are satisfied, i.e. the space 
of {qab, tt"'') such that H^{q, tt) = 0. We call this the constraint surface. Its dimensionality is 



24 



Gravitation quantique 



(12 — 4) • oo"^ = 8 • oo'^. The fact that the algebra of constraints is first class guarantees that the 
gauge transformations generated by the constraints preserve the constraint surface. We shall 
refer to the trajectories drawn by the gauge transformations as orbits. Points along one orbit 
correspond to the same physical configuration, only described in different coordinate systems. 
Hence, to select the physical degrees of freedom we have to divide by the gauge orbits, in 
a manner identical to what happens in gauge theories. Since the orbits span a dimension 4 
manifold at each space point, dividing by the orbits gives (8 — 4) • od'^ = 4 • oo'^. This is 
the physical phase space. It has four dimensions per space point, namely the theory has two 
physical degrees of freedom per space point (or simply two degrees of freedom, for brevity, 
with the space dependence tacitly implied), a result consistent with the linearized analysis. 
See e.g. [12] or [14] for more details. 



/ (9„6,^''^) = 12-oo3 / 


y- 

/ 
1 
/ 
1 
1 






__________ 




<^ 


» — — ' 


"^ {qab,^"^) =8- OO^ 
i? = iJ° = 


/ 


^ ( 


'~~~~~~~~~ 




f 


= 4-oo3 




/ 



This has far as the counting goes. However, in the case of the linearized analysis we are 
also able to identify the 2 degrees of freedom as the two helicities, and associate a physical 
trajectory to each point in phase space, thanks to the fact that we are able to solve the dyna- 
mics. Compare this to the mechanics of a point particle : the dynamics is given by (/'[^^.po] (^)' 
and the degree of freedom corresponds to varying the initial conditions (0O:Po); spanning in 
this way all possible physical evolutions. 

Therefore, if we want to know what the two physical degrees of freedom of general rela- 
tivity are, we need to control the general solution of the theory. This is a formidable task due 
to the high non-linearity of the equations, and in spite of the effort in this direction, still little 
is known. See [15] for a review of some attempts. 

There is a particular reformulation of the dynamics that plays a role in quantization. Let 
us go back to the Gauss-Codazzi equation (1.14). A simple manipulation gives the so-called 
first Gauss-Codazzi identity. 



n + K'' - Tt{KK) = 2ri^n''Gp 



(1.29) 



where G is the Einstein tensor. One can show (but we will not do it here) that the left hand 
side of (1.29) is equivalent to a linear combination of the constraints. This means that if g^^^, 
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satisfies the constraints H = iJ" = on any spacelike Cauchy hypersurface, than it also 
satisfy all ten Einstein equations G^^ = 0. This is a crucial point : the whole dynamical 
content of general relativity is in those four constraints. 

1.2.4 Dirac's quantization program 

The fact that the constraints include the whole dynamics is at the heart of the theory. 
It is heavily exploited by the approach to quantization proposed by Dirac, which is based 
on a definition of dynamical physical states as the ones annihilated by the constraints. The 
procedure can be schematically divided in three steps : 

(i) Find a representation of the phase space variables of the theory as operators in an auxi- 
liary "kinematical" Hilbert space Hkm, satisfying the standard commutations relations 

{•,-}^^ [•,•]; (1.30) 

in 

(ii) Promote the constraints to operators H^^ in Hkin ', 
(iii) Characterize the space of solutions of the constraints Hphys, 

H^'iP = Vt/. e -Hphys- (1.31) 

These steps should then be completed with an explicit knowledge of the scalar product in 
Hphys and a physical interpretation of the quantum observables. 

Dirac's procedure is more general than gravity : it applies to any fully constrained system. 
Let us now try to apply it to the ADM formulation of general relativity. We then look for a 
space of functionals carrying a representation of the quantum Poisson algebra 

[qabix),it'Hy)] = ^fiSl^b)SHx,y), (1.32) 

[qab{x),gcd{y)] = 0, 

[7t'^\x),7r^''{y)] = 0. 

Formally, we can proceed by analogy with better known cases, such as a scalar field theory, 
and consider a Schrodinger representation cjabix) = qab{x), tt°-''{x) — —ifi-g — ^rjy, acting on 
wave functionals 

^[qab{x)] (1.33) 

of the 3-metric. This procedure, which works well for the scalar field [4], encounters a number 
of difficulties when applied to the gravitational context. For instance, to define the (kinemati- 
cal) Hilbert space we need a scalar product, formally 

dgW]i^'[9] = {i^\i^')- (1-34) 

However, there is no Lebesgue measure on the space of metrics modulo diffeomorphisms 
that we can use to define dg. Without this, we can not even check that qab{x) and tt°'''{x) are 
hermitian, nor that qab{x) has positive definite spectrum, as needed to be a spacelike metric. 
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Let us ignore these and other similar issues, and try to proceed formally assuming that a 
well-defined Hkin exist. The next step is to promote the constraints (1.20) to operators, and 
characterize their space of solutions. Let us proceed in two halves, 

Tra rj TT n 

'Hfcm > T~(-Diff > Hphys- (1.35) 

Consider first the vector constraint. In the Schrodinger representation defined above, the 
smeared version gives 

H{Na)i^[qab\ = 2ih f d'^xVbNaP^ = 0, (1.36) 

after an integration by parts. This implies straightforwardly that 

■lp[qab + 2V(aiVb)] = ^Jj[qab], 

namely the solution of the vector constraint are those functionals of the metric invariants 
under diffeomorphism. This is very nice, as it realizes at the quantum level the correct action 
of the classical constraints. However, the space of solutions Hoiff is again ill-defined, since 
it inherits from the kinematical one the lack of measure theory or other means of control over 
it. 

For the Hamiltonian constraint we can write 



Hip[qab 



2 Vdctg oqab(x)dqcd(x) 



V'bah], (1.37) 



where the colon : means that an ordering of the operators needs to be prescribed. The situa- 
tion is more complicated for the Hamiltonian constraint, since it requires the definition of 
products of operators at the same point, notoriously very singular objects. The formal ex- 
pression (1.37) is known as the Wheeler-DeWitt equation. Even if manage to give a suitable 
ordering prescription and regularize the differential operator, the problem with the equation is 
that we do not have any characterization of the solutions, not even formally as for the diffeo- 
constraints above. And of course, again no clue on the knowledge of the physical Hilbert 
space and scalar product (see however minisuperspace models [16]). 

Loop quantum gravity is an approach to the problem which improves significantly the 
situation, and gives a number of answers to these open questions. The key to LQG and to 
such improvement is surprisingly simple ; instead of changing the gravitational theory or the 
quantization paradigm, we just use different variables to describe gravity. After all, we are 
familiar with the fact that not all choices of fundamentals variables work out as well when 
quantizing a classical theory. Consider for instance the harmonic oscillator. Classically, the 
most elegant description of the system is in terms of action-angle variables, which parame- 
trize the phase space as {(f), J} = 1 instead of {q,p} = 1. However, it is more convenient to 
quantize the system using the {q,p) variables than the action-angle ones, which require extra 
care in constructing the operator algebra and dealing with unitarity. 

We now introduce the variables that allows us to reformulate general relativity in a way 
more amenable to Dirac's quantization procedure. 
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1.3 Appendix 

In this Appendix we report some calculations which are used in the main text. 

Proof of the Gauss-Codazzi equation. We start the computation from the characterization of 7?. in terms 
of covariant derivatives V defined as the projection of V on E via q^ 

2Vi^V,]Wp = -■wp,n''p^, (1.38) 

V^V.^p = V„ [q:'q'p'V,,w,,^ = q^'qfqfV^, {q:'"q;„^ .'Wp,^ (1.39) 

it's useful to notice that 

q^'q^'y^^'Qt' = Qm'?^' V^' {gt' + n'n,,) = (1.40) 

= <'?:' V^,5^, + g,'l'gr'n,, V^m" + q'^' q^' n''V^,n,, = (1.41) 

= "'''?M 9^ V^/n^/ = K^.^n'' (1.42) 

we used the definition of extrinsic curvature, orthogonality (q'^ n^i = 0) and metric compatibility 

(V^-fl^, = 0). 

V^V^Wp = g,*^ ql ql V^/ V^-Wp- + g^ n A'pi.V^'iUp/ + q^ K^pii'' V^>Wp> (1.43) 
in particular if we consider u; G E thus n'' Wpi = 

nf V^iWpi — V^in'' Wpi — WpiV^in'' — —WpiV^in^ (1.44) 
Antisymmetrizing in a,b 

V[^V^]Wp = q^^q^]qp Vp'V^'Wp, + q^ n" K[pr,]Vr,'Wp> - q[„ K^]pWp,\/^,n'' = (1.45) 

= 9^ 9^ 1p V[^'V^/]w;p/ - ii'[''^ii'p]pU)p/ = (1.46) 

= -^q^' q,' qi ^p'p, ,, Wp,, ~ k(^ K^^pWp, (i .47) 

or in terms of the Riemann tensor using (1.38) 

7^?M.»P' = [€iUi'Kn^,,. + 2;ff>p,p) Wp, (i.48) 

this equation state for w G E for make it true for all « G 7V/( it is sufficient to take Wp — q^ Vpi 

7?-pM^ = 9m 1^ 1p q"cy'K'p'v' + S-K'f^ii'Mip (1.49) 

Or equivalently 

Tlappu = q*^ ql q^ qZ R^'p'p'u' + "iK^y^Kp^p (1.50) 

7^ = q'^'q^'q^p ql qp ql R^'p'^'v' + q^^q^^K^^Rpp - q'"^q''''KapK^p = (l.5l) 

= q^'^'q^'^'R^'p^p,,, + Tr(A-A-) -K^ = (1.52) 

= i? + 2n"'n'''i?pv/ + n'''n'''n'''n''' R„,p,p,^, + Tr{KK) - K'^ (1.53) 
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the term nf' n'^ n" n'' Ra'p'y.',y vanishes because of the symmetries of the Riemann tensor, using 



the fact that n"^ n'' gpi^i = — 1 we can write 



7^ = 2n" n" ( Rp,,y - ^^R) + Tr{KK) - K^ (1.54) 



That is the first Gauss-Codazzi equation. 



Constraint Algebra. We define the Poisson bracket as 

We start the computation from the diffeomorphism constraint 

H{N) = f d^xN'Hi (1.56) 

we want to show that it generates infinitesimal spatial diffeomorphism (it's sufficient to do it on the 
phase space variables) 

[qativ), H{N)] = J d\N' ^^ = 2V(„iV6) = L^qat (1.57) 

in fact from the definition of the Lie derivatives 



Note that the functional derivative is easy to compute with a by part integration. The last term vanishes 
for metric compatibility and in the Lie derivative you can substitute d with V. 

I- ' J Oqabiy) ^ \y^ 

V^ V^ ^/9 

Note that tt"'' is a tensorial density so its lie derivative has to be defined as above. Given A an arbitrary 
function of q and tt using Leibniz we can compute the Poisson bracket easily 



{a,h{n)} 



5A SH{N) SH{N) SA 



i^{,..,«,N,}-p,«,N)}|l 



Sq, 

SA _ ^b SA 

Sqab Sir 



Cf^qab — C^TT — — ^ — -CjyA 
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In particular for H and Ha as function of h and n we deduce 



.Ik + Ik 



^^Hb = y/qCf}—- + —-Cf^^ ■- 



= N^daHb + N^Hty/^da-^ + HadbN" + HbVaN" = 

= N^daHb - N^Htdai^)-^ + HadtN" + HtVaN" 

= N^daHt + HadtN^ + HtiVaN'^ - ^daVq) = 

= N^daHt + HadbN" + HtdaN" 



In fact 



VaN" =daN''+rtaN'' 

M"- 1 

—pda^fq = T:N"-q''^daqab = KaN'' 

Where the last equality is due to metric compatibility. We can conclude that 

N^Ha, hA^ -N^daHb - HadtN'' - HtdaN" (1.58) 



That is the smeared version of 

{Ha{x), Ht{y)} = Ha{y)dbSix -y)- Ht{x)d'j{x - y) (1.59) 

Moreover also the second Poisson bracket is easy to compute using the fact that H{N) is the generator 
of infinitesimal diffeomorphism along N. 



H_ H_ 



Cj^H = ^r^— + -=Cf,^ = (1.60) 



H H 1 . 



= yfqN^da—p + —p{-y/qq"^Cf}qad) = 



= N''daH + ^N^HdaA- + ILh^q-\^^y^N^ = 

= N-'daH + HdaN"" + ^N^HA^q'^daqcd + HVI^N'' = 

2 V? 

=0 

putting all together 

'' N^Ha, h\ = -N^daH - HdaN" (1.61) 



that is equivalent to 

{Ha{x),H{y)} = Hda&{x-y) (1.62) 
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To ultimate the computation of the second Poisson bracket we need some preliminary results. Given 
an arbitrary function /(vr, h) 

{{f,H{m)} ,H{N2)] = C^^^C^^J (1.63) 

using the Jacobi identity 

{/, [H{N^),HiN2)}} = [[H{N2),f},H{Ni)} + {[f,H{Ni)},H(N2)} = 

If we smear the (1.61) we obtain : 

' ' N^Ha, f Nh\= f CffNH (1.64) 



(1.65) 



I f N^Ha, f NSHa\ = f [Nl,N2r Ha 

Finally we compute the last Poisson bracket, if N ^ N' (if they are equal it vanish) 



{HjNlHjN')} = /d-. f f , f If ; - f (f j f y) (1.66) 

We notice that all the algebraical terms in q and tt simplify each other. In fact 

SH[x 



Sqab{y) 
SH(x) 



S-^abiy) 



= t{x)ab5{x-y) (1.67) 

algebraic 
part 



= U{xrS{x-y) (1.68) 

algebraic 
part 



{H[N)\algebraic , H [N )\algehraic \ 
part part J 



d'xNf'^bN'ff - N'ffNf^b = (1.69) 



The non algebraic part in q are in R conversely the Hamiltonian constraint is algebraic in vr. 







non non 

t«^ivj,ii^iv Ji yax ^^^^^^^ ^^^^^^^ ^^^^^^^ ^^^^^^^ 



We compute (ignoring algebraic terms for simplicity) 
Sg f NH = - f N^q'^'SRat = 



= - / Ny:^ [y aV bSqcdq"" q"^ - q'"^V''VbSqcd) = 



(1.70) 
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Substituting this variation in the Poisson bracket the first term becomes 
SJNHSJN'H 



Sqab Sn"'' 



y^ (VV'iV - V'(iVg"'')) 2N' ( ^ 



T^Qab 



V^ ^/g(3 - 1) 



= 2 / (V"V''iV)7r„6iV' + 



3-1 3-1 



3- 1 

^ 11 ttN'v'^N 



3-1 



2 / yg(V"V'iV)^iV' = 2 A ygV^ATV 



^/9 



(1.71) 



summing also the second term 
[1.70] = 2 / Vg 
= 2 



V^WV" ( ^^^ - V^iY'V ^^""^ 



^/5 



^ 



(1.72) 



i L \/g V\/g/ \/g 

-Af'V'AT + iVV^AT') 2^V" ('^'j = f (n\/''N' - iV'V'iv) ffa 



lb at'^yjo. I '^ab 



V^ 



2 Tetrad formulation 

A tetrad is a quadruple of 1-forms, e'Ax), 7 = 0,1,2,3 such that 



(2.1) 



By its definition, it provides a local isomorphism between a general reference frame and an 
inertial one, characterized by the flat metric ?//j. A local inertial frame is defined up to a 
Lorentz transformation, and in fact notice that the definition is invariant under 



eiix) ^ ?/(x) ^ A'jix)e'(x). 



(2.2) 



This means that the "internal" index / carries a representation of the Lorentz group. Contrac- 
ting vectors and tensors in spacetime with the tetrad, we get objects that transform under the 
Lorentz group, e.g. ef^n'' = n^ . The tetrad thus provides an isomorphism between the tangent 
bundle of 7W,r(7V/() = Up (P>rp(-^))' and a Lorentz principal bundle F = {M,S0{3,1)). 
On this bundle we have a connection cul'' , that is a 1-form with values in the Lorentz algebra, 
which we can use to define covariant differentiation of the fibres. 



D^v'{x)^df,v'{x)+^ij{xy{x). 



(2.3) 



This is the analogue of the covariant derivative V^ = 9^ + F^ for vectors in T{Ai). We can 
also define the derivative for objects with both indices, such as the tetrad. 



T^^ei 



dA+^ij4~^iA- 



(2.4) 
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As the Levi-Civita connection T{g) is metric-compatible, i.e. V ^^gup = 0, we require lj^ to 



be tetrad-compatible, i.e. P^e^ = 0, and call it spin connection. This implies 

diA) + ^Ij<) = K,/p^ ^^A + <■< = KAp ^ 0' (2.5) 

where we separed the spacetime indices into their symmetric and antisymmetric combina- 
tions, and used the fact that the Levi-Civita connection V{g) has no antisymmetric part. 

From these equations we immediately obtain the following relation between the spin and 
Levi-Civita connections, 

^Ij = 4v^e:;, (2.6) 

as well as the fact that the spin connection satisfies 

d^e^ = de^ + u^j he-' ^ (a^e^ + uljci) dx'^ A dx"" = 0. (2.7) 

This equation is known as Cartan's first structure equation. Here we introduced the exterior 
calculus of forms, with d the exterior derivative, d^^ the covariant exterior derivative, and A 
the wedge product. See [17] for an introduction to this formalism. 
Given the connection, we define its curvature 

F" = dw" + uj^K A Lo^^, (2.8) 

whose components are 

Using the solution uj(e) given in (2.6), an explicit calculation gives 

F;^;^ {cu (e)) = e'Pe''-R^,p4e), (2.10) 

where Rp,vpa{s.) is the Riemann tensor constructed out of (the metric defined by) the tetrad 
efj. This relation is known as Cartan second structure equation. It shows that general relativity 
is a gauge theory whose local gauge group is the Lorentz group, and the Riemann tensor is 
nothing but the field-strength of the spin connection. 

Cartan second structure equation. Starting from the definition of F in (2.9) and inserting (2.6), we have 

+ epLJ^e^e,^ d^e + epd^ej^e^ T^^e + eprj^e^e^ r^i,e ' - (^ o !/). 

Next, we use e^e^ = 5]^ and Cpd^e'^ — —dp{ep)e''i^ to rewrite this expression as 

Fpi^ = Cpe'-'dp {V^,,) + e'pc'^-'r'^pr^^ + dmue'pd^e'"'' - dpe'pd^e''--' + 

+ dpCpV^^e"'' - dpCpV^^e"'' + epV^^dpc"'' + CpT'^^pdve"'' - {^ -k^ v) 

Relation between the determinants. The determinant g of g^v is related to the determinant of the tetrad 
e by the simple relation 



g = -e\ (2.11) 
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This expression can be easily derived recalling Caley's formula for the determinant of a matrix, 

g = detp^„ = ^e'"''"^e"''^''pM<:«5^/55p7fl<T«- (2.12) 

If we substitute the expression of g^„ in terms of tetrads we get 

J , / / J \ 1 uuprr aHj6 I J K L M N OP 

g = det [e^e^riij j = — e e e^earjije^ epriKLep e^ T)MNe„ eg -gop = 

1 nvntr I K M O aR-iS J L N P 

= —£ e^e^ Bp e^e e^e^e^ eg riij-gKLVMNrjop = 

1 2 IKMO JLNP 2 

— j^e e £ gijr]KLgMNrjOP = — e 



2.1 The action in terms of tetrads 

The Einstein-Hilbert action can be rewritten as a functional of the tetrad in the following 
way (recall we take units IGttG = 1), 

SEH{el) = \jijKL I e' he^ A F^^ {u{e)) . (2.13) 

On top of the invariance under diffeomorphism, this reformulation of the theory possesses an 
additional gauge symmetry under local Lorentz transformations. 

Rewriting the actions in terms of tetrads. Explicitly using (2.10) and (2.11) 

SEH{gp.u{e)) = / d'^x^f^g^'^R^i, = / A'^xee^^e"' Rf_,p^cre''je'''' = 

= / d*xee''je''jFlp {Lo{e)) = I A^x-eukle''''"''^ e'^e^^F^p {u[e)) 
= f^eijKLe Ae-' AF'^'^iuie)) 

A fact which plays an important role in the following is that we can lift the connection to be 
an independent variable, and consider the new action 

S{ei,u;'/) = \eijKLJ e' A e' A F^^ H . (2.14) 

Although it depends on extra fields, this action remarkably gives the same equations of motion 
as the Einstein-Hilbert one (2.13). This happens because the extra field equations coming 
from varying the action with respect to w do not add anything new : they simply impose the 
form (2.6) of the spin connection, and general relativity is thus recovered. 

Deriving the fields equations. One first verifies the Palatini identity 

5„F^^H =d„5a;'^''. (2.15) 
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Then, the variation gives 

SuS = ^eijKL f e' Ae-' A d„5a;^^ = -^Sijkl f d^ (e' A e'') A 5a;^^ (2.16) 

after an integration by part. Imposing the vanishing of the variation, we obtain the field equation 

euKLc' Ad^e-' = 0. (2.17) 

If the tetrad is invertible this equation implies duie'^ — 0, which in turns implies (2.6) in terms of the 
Levi-Civita connection of the metric associated with e^. 

As it gives the same field equations, (2.15) can be used as the action of general relativity. 
Notice that only first derivatives appears, thus it provides a first order formulation of gene- 
ral relativity. Furthermore, the action is polynomial in the fields, a desiderable property for 
quantization. On the other hand, there are two non-trivial aspects to take into account : 

- The equivalence with general relativity holds only if the tetrad is non-degenerate, i.e. 
invertible. On the other hand, (2.15) is also defined for degenerate tetrads, since inverse 
tetrads never appear Compare the situation with the Einstein-Hilbert action, where the 
inverse metric appears explicitly. Hence, the use of (2.15) leads naturally to an exten- 
sion of general relativity where a sector with degenerate tetrads, and thus degenerate 
metrics, exists. 

- If we insist on the connection being an independent variable, there exists a second term 
that we can add to the Lagrangian that is compatible with all the symmetries and has 
mass dimension 4 : 

^/jifLe^Ae-^AF^^H, (2.18) 

where Sjjkl = Si[kSl]j- This term is not present in the ordinary second order metric, 
since when (2.6) holds, 

SijKLc' A e'' A F^^(c^(e)) - e^'^'i^^.p^Ce) = 0. (2.19) 

Adding this second term to (2.15) with a coupling constant I/7 leads to the so-called Hoist 
action [18] 

S (e, uj) = (^eijKL + -Sijkl) f e' A e' A F'^"^ {uj) . (2.20) 

Assuming non-degenerate tetrads, this action leads to the same field equations of general 
relativity, 

^" ^ei^t^e''", Gp,(e)=0. (2.21) 

This result is completely independent of the value of 7, which is thus a parameter irrelevant in 
classical vacuum general relativity. It will however turn out to play a key role in the quantum 
theory, where it is known as the Immirzi parameter. ^ 



Deriving the fields equations for the Hoist action. The vanishing of the variations give 



-eijKL + -SijKL I e' A dt^e'^ = 0, 
2 1 J 

\eijKL + -SuKi\ e-' A F^^ {uj) = 0. 
2 1 ) 



(I.IT) 
(2.23) 



2. The Immirzi parameter becomes relevant also at the classical level if source of torsion are present 119, 20]. 
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For invertible tetrads, the first one is again uniquely solved by (2.6). Substituting this solution into the 
second equation we get 

e G„ H — oijkle eu-TaS = e G^ H — oijkle e„e e '^ HaSvp = 

= e G^ + — e*^ ( e d^ - e '^d^j iJ^i^p = e Gy e ^ e R^s^^ = e G,^ = 0, 

which is equivalent to G^u, = thanks to the non-degeneracy of the tetrad. In the last step we used the 
first Bianchi identity e'^'^^^-Rai^,. = 0. 



2.2 Hamiltonian analysis of tetrad formulation 

For the Hamiltonian formulation we proceed as before, assuming a 3 + 1 splitting of the 
space-time {M. = R x S) and coordinates (t, x). We introduce the lapse function and the 
shift vector (N, N"-) and the ADM decomposition of the metric (1.11). It is easy to see that a 
tetrad for the ADM metric is given by 

ei = e^r^ = Nn' + TV^e^, 5,,e>^, = gab, * = 1, 2, 3. (2.24) 

The "triad" e^ is the spatial part of the tetrad. As before, we want to identify canonically 
conjugated variables and perform the Legendre transform, but we now have two new features 
which complicate the analysis. The first one is the tetrad formulation, which in particular has 
introduced a new symmetry in the action : the invariance under local Lorentz transformations. 
As a consequence, we expect more constraints to appear, corresponding to the generators 
of the new local symmetry. The second one is the use of the tetrad and the connection as 
independent fields. Therefore, the conjugate variables are now functions of both e^ and lo^^ 
(and their time derivatives), as opposed to be functions of the metric gab only. 

The consequence of these novelties is a much more complicated structure than in the 
metric case. In particular, the constraint algebra is second class. However, there is a particular 
choice of variables which simplifies the analysis, making it possible to implement a part of 
the constraint and reducing the remaining ones to first class again. These are the famous 
Ashtekar variables, which we now introduce. ^ 

To simplify the discussion, it is customary to work in the "time gauge" e^ n^ ~ 5q, where 



el = {N,Q)^el^{N,N'^el). (2.25) 

The crucial change of variables is the following : we define the densitized triad 



Et - ee,« = -e,,ke'''''elet (2.26) 



and the Ashtekar-Barbero connection 



^:.-7c^a" + ^e>f. (2.27) 



3. On the general analysis with the second class constraints see [21]. 
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These variables turns out to be conjugated. In fact, we can rewrite the action (2.20) in terms 
of the new variables as [21, 13] 



S{A,E,N,N'') 



1 



dt I d^x 

s 



A^E^ ~ A^DaE^ - NH - N^Ha 



where 



Gj = DaEf = daE^ + EjuAiE''' 



^^ - -K^^) 



H = 



FL - il' 



7 



K^Gi, 



j7nnJ-^a 



detE 



1+^2 ^a 

—-^G'da ' 



detE' 



(2.28) 

(2.29) 
(2.30) 

(2.31) 



The resulting action is similar to (1.17), with {A, E) as canonically conjugated variables, as 
opposed to (g, tt). Lapse and shift are still Lagrange multipliers, and consistently we still 
refer to H{A, E) and Ha{A, E) as the Hamiltonian and space-diffeomorphism constraints. 

The algebra is still first class. The new formulation in terms of tetrads has introduced the 
extra constraint (2.29). The reader familiar with gauge theories will recognize it as the Gauss 
constraint. Just as the H^^ constraints generate diffeomorphisms, the Gauss constraint gene- 
rates gauge transformations. It is in fact easy to check that E^, and A]^ transform respectively 
as an SU(2) vector and as an SU(2) connection under this transformation. 

The algebra generated by the Gauss constraint. We define the smearing of the Gauss constraint 

G[A)= I A^xGr{x)A\x). (2.32) 

Its Poisson bracket with the canonical variables give 

d:'xk'{x)G,{x),E^{y)^ = j d\k'{x) [c 



{/■ 



\d,E° + em,^ATE''',E^{y)] 



Hence, 



d'xK'{x)G,{x),A\{y) \ = / d'xK^x) [d^E] + Sm^r^AT E""^ , A 
= 7a.A'(j/) + 7£.n,.A^(j/)^;"(y). 



'ai.y)} 



+ 



[G{k),daEny) + e,.uA'aE''\y)] = 

= ier,„da (A^(y)i5""(j/)) -7£..feA^a(y)£MnA'(y)£°"(j/)- 

+ 7e,.fei?'"' (da^'{y) + e^t,K\y)A":{y)) = 

= 7eij„A^ {y)daE"\y) + 7 {eijkSe„k - einkeijk) A\y)E'"'(y)Ai(y) = 

= jeukA'iy)daE''''{y)+-,eukA\y)e,knE''"iy)Ai{y) = ^sukA'{y)Gkiy). 
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Finally, smearing also tlie second term in the above bracket we get 

{G(Ai),G(A2)} = |g([Ai,A2]), (2.33) 

wliicli we recognize as tlie su{2) algebra structure equations. 

We should not be surprised of the appearance of this extra constraint. When we use the 
tetrad formalism, we introduce a new symmetry in the theory, the invariance under local 
gauge transformation. The Gauss constraint is there to enforce this invariance at the canonical 
level. On the other hand, it might be more puzzling that although the local gauge invariance 
of the covariant action was the full Lorentz group, the Legendre transform (2.28) is only 
invariant under SU(2). 

The origin of this puzzle Ues precisely in the change of variables (2.26-2.27). The Ashtekar- 
Barbero connection is an SU{2) connection, not a Lorentz connection. In particular, it is not 
the pull-back of the space-time Lorentz connection [22]. It should then be seen as an auxiliary 
variable, useful to recast the algebra in a first class form. 

Summarizing, in this formulation of General Relativity the theory is described by an 
extended phase space of dimension 18 • oo"^ with the fundamental Poisson bracket 

{A^,{x),E]{y)]=^5'j]5\x,y) (2.34) 

This new internal index i corresponds to the adjoint representation of SU(2), and we can 
recover the old 12 • oo"^ dimensional phase space on the constraint surface d — Q dividing by 
gauge orbits generated by G. This group SU(2) should be seen as an auxiliary local symmetry 
group, since the connection with the original Lorentz group of the tetrad formulation is hidden 
in the change of variables (2.26-2.27). It is only for the special case 7 = 1 (which is the 
original one introduced by Ashtekar) that the relation is manifest ; in that case, the SU(2) 
corresponds to the self-dual subgroup of the Lorentz group. Notice that this case also leads to 
a simplification of the Hamiltonian constraint (2.31). On the other hand, the variables are now 
complex, and to recover general relativity one needs to impose reality conditions. These are 
particularly difficult to deal with at the quantum level, and for this reason most developments 
in LQG have focused on 7 real. This is what we do also in this review. 

2.3 Smearing of the algebra 

The next step is to smear the algebra (2.34), as we did previously with the ADM variables. 
This is needed in order to proceed with the quantization. At this stage, the different tensorial 
nature of Aa and E"- plays a key role. If we were doing gauge theory on flat spacetime, we 
would not pay attention to the different indices, and just smear the connection and the electric 
field with the same type of test functions. On the contrary, now that the non-trivial geometry 
of spacetime is our main goal, we have to be careful. In fact, a brief look at (2.26) shows that 
the densitised triad is a 2-form. Hence, it is natural to smear it on a surface, 

E,{S) = I UaE^d^a, (2.35) 

Js 

where Ua = Sabc^-^- is the normal to the surface. The quantity Ei{S) is the^Mx of E 
across S. 
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The connection on the other hand is a 1-form, so it is natural to smear it along a 1- 
dimensional path. Recall that a connection defines a notion of parallel transport of the fiber 
over the base manifold. Consider a path 7 and a parametrization of it x°'{s) : [0, 1] -^ E. 
Given a connection A^ we can associate to it an element of SU{2) Aa = A^^Ti where r^ are 
the generator of SU{2). 

7 / Then we can integrate Aa along 7 as a line integral. 




A= / dsAl{x{s))^T^T,. 
7 "'o 



ds 



Next, we define the holonomy of A along 7 to be 



/i^ = P exp ( / Aj , 
where V stands for the path-ordered product. That is, 

h^^Yl /// A{j{si))---A{j{sn))dsi---dsn 

"="l>s„>--->si>0 

where we parametrized the line with s E [0, 1]. 



(2.36) 



(2.37) 



(2.38) 



On the definition of the holonomy. More precisely, we call holonomy the solution of the differential 
equation 

^/i^(i)-/i^(i)A(7(t)) = 0, ft^(0) = l. (2.39) 



If we integrate the equation by iteration we have 

hj(t) = l+ A('y(s))h^{s)ds = 



(2.40) 



= 1+/ A{j{si))h^{si)dsi + / A(7(si))A(7(s2))/i^(s2)dsids2 

Jo Jo Jsi 



therefore formally 



h-,{t) = Yl /// A{j(si)) ■ ■ ■ A('y{s„))dsi ■ ■ ■ dsn 

"=9>si>--->s„>0 



(2.41) 



To complete the proof, one needs to show that the series is well defined. Indeed, it converges with 
respect to sup norm 



11^711 <E 



t" 



|A(7(si)) • . ■ A(7(s„))|| dsi . . . ds„ < V Ij All" '-. (2.42) 

-^ — ^ rji 



"t>si>--->s„>0 



For further reference, let us also notice that the terms of the series can be written as integrals over square 
domains (si, . . . , s„) £ [0, t]", instead triangle domains f > si > ■ • ■ > s„ > 0. This gives 



h-y(t) = Vexp 



A 



UIIL 



P(A(7(.i))...A(7(s„)). 



(2.43) 
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Let us list some useful properties of the holonomy. 

- The holonomy of the composition of two paths is the product of the holonomies of 
each path, 

hpa = hpha. (2.44) 

- Under a local gauge transformations g{x) gSU(2), the holonomy transforms as 

^^ 9s(^)h^9t(^Y (2.45) 

where 5(7) and ^(7) are respectively the source and target points of the line 7. 

- Under the action of diffeomorphism, the holonomy transforms as 

h^ {dp* A) = h^o^ {A) . (2.46) 

- The functional derivative with respect to the connection gives 



5h^[A] 



ii°(5'^-'(7(s),a;) Tih^ if 2; is the source of 7 

= ■( ii°5'^^)(7(s),x) /i^Tj if X is the target of 7 (2.47) 



"■^ ' y x°-5^^>{^{s),x)h-y{0,s)Tih^{s,l) if X is inside 7 



Basic proofs. If not otherwise specified, we assume for simplicity that the source is to = 0. For two 
composable paths a and /?, we define the composition 

^""-X ait-S), in€[S,T + S] ^^-^^^ 



If we split the integrals as J = J + J„ , (2.43) reads 



s 

S rS nT nT 



^\.] I dti--- dti I dti+i ■■■ dt„. (2.49) 



° ,=() VJ ■'0 Jo Jo Jo 

This allows us to split the path ordering as follows, 

P {A(/3a{h)) • • • A{l3a{t„))) = P (yl(/3(i,+ i)) ■ ' ■ A{l3{t„))) P {A{a{h)) • • • A{a{U))) . (2.50) 
Hence, the holonomy along the composite path is 

^/^"-EE -.f -V / V {Am.+^)) ■ ■ ■ Am„))) PiAia{t,))---Aia{U))) 

= E 4 / ^ (Mm)) ■ ■ ■ Am,))) t^ f P {A{a{t,)) ■ ■ ■ A{aiU))) 
i^o g! Jam *! JniS) 

= hpha, 

which proves (2.44). In the second step, we defined q = n — i. 

To prove (2.45), let us introduce families of vectors u{t) and w{t) such that 

w{t) = u(i)g(7(t)) = h^{t)g(jit))u(0). (2.51) 
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w{t) satisfies the following differential equation, 

—w{t) = hgu{0) + hgu{0) = hgg-^gu{0) + hgg-"^ Agu{Q) 

= hg {g-'g + g'^Ag) u{Q) = {gg-' + g-^Ag) ghu{0) = A'w{t), 

which implies that w{t) = h^t)w(0) = h?,{t)g{'y{0))u{0). Comparing with (2.51) we find (2.45). 
The action of a diffeomorphism on the line integral of the connection is given by 

^ <t>A^ j A^{ci>o^{t))d.^^{-i{t))Y{t)dt = j A^{^o'y{t))j^{ci>o^Y{t)dt, (2.52) 

which implies immediately 

h^ {(t>A) = h^a-, (A) . (2.53) 

Finally, to prove (2.47), let us compute the differential equation satisfied by ,"! ■ ,"' in the two 
cases of x being inside the path or at one boundary. If x is inside 7, from (2.39) we get 

d 5hj(to,t) Sh.j{to,t) 



dt SAl{x) 5Al{x) 



-Ai-yit)) = 0, (2.54) 



which is solved by ^^ ■ ?\ = h{to, s)Tix"'h{s, t). Then, by the Leibniz rule we also know that 



SAl{x) 



T#V = %¥^^^.(M) + ^.(to,.)%|f4 = Mto,.)r.i"M.,t). (2.55) 

5Al{x) oAl,{x) oAl[x) 



From this we can argue that if x is at one boundary we must have 

5hj{to,t) _ J ^x'^Tih^{to,t) if 2; is the source of 7 



5Ai{x) \ 2^°'hi{to,t)Ti if s is the target of 7 



(2.56) 



2.4 Summary 

In this section we have taken the two key steps needed to prepare general relativity for the 
loop quantization. The first step was to reformulate the theory in terms of the tetrad field and 
an SU(2) independent connection, the Ashtekar-Barbero connection (2.27). The second step 
was to regularize the resulting Poisson algebra using paths and surfaces, instead of the all of 
space as in traditional smearings (cf. what was done in the ADM formulation). The resulting 
smeared algebra of h^[A\ and Ei{S) is called holonomy-flux algebra. It provides the most 
natural regular (i.e. no delta functions appear) version of the Poisson algebra (2.34). 

3 Loop quantum gravity : Kinematics 

The formulation of General Relativity in terms of the Ashtekar-Barbero connection and 
the densitized triad let us to talk about General Relativity in the language of a SU{2) gauge 
theory, with Poisson brackets (2.34) and the three sets of constraints 

Gi = Gauss law 

Ha = Spatial diffeomorphism invariance 

iJ = Hamiltonian constraint 
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The difference with a gauge theory is of course in the dynamics : In gauge theory, after 
imposing the Gauss law, we have a physical Hamiltonian. Here instead we still have a fully 
constrained system. 

Let us first briefly review some basic steps in the quantization of gauge theories. In a 
nutshell, the usual procedure goes along these lines : 

- Use the Minkowski metric to define a Gaussian measure 5 A on the space of connections 
modulo gauge-transformations. 

- Consider the Hilbert space L2 (A, 5 A) 3 ^^J[A] and define the Schrodinger representa- 
tion 

Ali,[A] = A>[A], (3.1a) 

E^i^lA] ^ ~thjjj-^p[A], (3.1b) 

which satisfies the canonical commutation relation, 

'Al{x),&^iy)] =ih^5l5]5\x,y). (3.2) 

- Impose the Gauss law constraint, GiijjlA] = 0, which selects the gauge-invariant states. 

- Study the dynamics with the physical Hamiltonian, 

idt^lj[A] = H^[A]. (3.3) 

In practice, the easiest way to deal with the dynamics is to decompose A in plane waves 
and the Hilbert space into a Fock space, L2{A,5A) = T-Lpock = ®„^n> the direct 
sum over n-particles spaces. 

The key difference in the case of general relativity is that we do not have a background 
metric at disposal to define the integration measure, since now the metric is a fully dynamical 
quantity. Hence, we need to define a measure on the space of connections without relying on 
any fixed background metric. The key to do this is the notion of cylindrical functions, which 
we introduce next. 

Commutator. 

\K{x),E]{y)\ ,p[A] = -ih^^^[x)-^^p[A]+ih^-^^^{x)^P[A] 

= ^^^VTfM'^t^] = ih-y5lS]S\x,ymA] 
SAiiy) 



3.0.1 Cylindrical functions and the kinematical Hilbert space 

Roughly speaking, a cylindrical function is a functional of a field that depends only on 
some subset of components of the field itself. In the case at hand, the field is the connection, 
and the cylindrical functions are functionals that depend on the connection only through the 
holonomies he [A] = V exp ( J^ A) along some finite set of paths e. Consider a graph F, 
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defined as a collection of oriented paths e C S (we will call these paths the links of the 
graph) meeting at most at their endpoints. Given a graph F C S we denote by L the total 
number of links that it contains. A cylindrical function is a couple (F, /) of a graph and a a 
smooth function / : SU{2)^ — !• C, and it is given by a functional of the connection defined 

as 

{A I F,/) = ^(rj)[A] = fiKAA],.. .,K,[A]) e Cyh (3.4) 



where e^ with 4 = 1, 



TJ)=4'irj)[A]=f{KMl---,heAA])eCylr 
. . ,L are the links of the corresponding graph F. 




Collection of paths F = {ei 



This space of functionals can be turned into an Hilbert space if we equip it with a scalar 
product. The switch from the connection to the holonomy is crucial in this respect, because 
the holonomy is an element of SU(2), and the integration over SU(2) is well-defined. In parti- 
cular, there is a unique gauge-invariant and normalized measure d/i, called the Haar measure. 
Using L copies of the Haar measure, we define on Cylr the following scalar product, 

{^Pirj)\^(rj'))= fUdhJiKM], ■ ■ ■,h,,[A])f\h,,[A], ■ ■ ■ ,/ie,[^]). (3.5) 

-' e 

This turns Cylr into an Hilbert space "Hr associated to a given graph F. 

Next, we define the Hilbert space of all cylindrical functions for all graphs as the direct 
sum of Hilbert spaces on a given graph. 



n. 



m 



res 



(3.6) 



The scalar product on Hkin is easily induced from (3.5) in the following manner : if Tp and 
ip' share the same graph, then (3.5) immediately applies. If they have different graphs, say Fi 
and F2, we consider a further graph F3 = Fi U F2, we extend /i and /2 trivially on F3, and 
define the scalar product as (3.5) on F3 : 



(V'(ri,/i) I ^{r2,f2)) = (Viur^./i) I ^(riur^./^)) 



(3.7) 



The key result, due to Ashtekar and Lewandowski, is that (3.6) defines an Hilbert space over 
(suitably generalized, see [25] for details) gauge connections A on S, i.e. 



Ukin = L2[A,dHAL]- 



(3.8) 
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The integration measure dfiAL over the space of connections is called the Ashtekar-Lewandowski 
measure. What (3.8) means is that (3.7) can be seen as a scalar product between cylindrical 
functionals of the connection with respect to the Ashtekar-Lewandowski measure : 



d/iAL V'(ri,/i)(^)V'(r2j2)(^) = (V'(ri,/i) I i^ir,.h)) ■ (3.9) 

Now that we have a candidate kinematical Hilbert space which does not require a back- 
ground metric, let us look for a representation of the holonomy-flux algebra on it. To that 
end, it is convenient to introduce an orthogonal basis in the space. This can be done easily 
thanks to the Peter- Weyl theorem. It states that a basis on the Hilbert space L2{G, d^Haar) 
of functions on a compact group G is given by the matrix elements of the unitary irreducible 
representation of the group. For the case of SU{2), 

fig) = E /™nA<»n(5) i^y ' ^' ■ ■;■ (3.10) 

"^ J 1 • • • 1 J 

1 

where the Wigner matrices Dmn{g) give the spin-j irreducible matrix representation of the 
group element g. This immediately applies to Hr, since the latter is just a tensor product of 
L2{SU{2), AjiHaar)- That is, the basis elements are 

{A I r;je,me,ne) = D\i\i^{K,)...D^^-XS^-^)^ (^-H) 

and a function V'(r,/) [^] G ^r can be decomposed as 

i^{Tj)[A]= J2 /^"::::frn„,,u,...,„„^HiL(/^e,[^])...i?H:L(/^e„[A]). (3.12) 

On this basis, we can give a Schrodinger representation like (3.1) for the regularized 
holonomy-flux version of the algebra. Consider for simplicity the fundamental representation, 
he = D^^^{he\ The holonomy acts by multiplication, 

h^[A]he[A] = h^[A]he[A\, (3.13a) 

and the flux through the derivative (2.47), 

E,[S)K[A] - ~ih^ f d^^na J.l'ff.. = ±ihjh,,[A]nh,,[A]. (3.13b) 

Here ei and 62 are the two new edges defined by the point at which the triad acts and the sign 
depends on the relative orientation of e and S. The action vanishes, E[S]he-^ [A] = 0, when e 
is tangential to S" or e n 5* = 0. 

Action of the Flux. In the case e U 5 = p recalling (2.47) we have 

Ei{S)K =-ih-y f dVn, / dsx''{s)S'\p,x{a))h,,rih,^ = (3.14) 

J s Jo 

- ih-yi-iy^P^h^iO, s)nhe[s, 1) (3.15) 
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where (— l)°(f ) = ±1 depending if ai, (T2, s is a right or left handed coordinates system 

1 i^t" r)r'' f)r'^ f 

daidaadse^b,-— s\x{s),x{a)) = ± / A\S{x) = ±1 (3.16) 

OCTl d(J2 OS J 

If there are no intersection the integration of the 5'^ (p, x{a)) is 0. In the case of e tangent to S we can 
compute the action of i? as a limit of a double intersection. Because the two contributes have relative 
opposite sign the limit is trivially 0. 

Consider now the action of the scalar product of two fluxes acting inside the link, 

E,{S)E\S)K[A] = -h^l^ h,,[Aynh,,[A]. (3.17) 

On the right hand side, we see the appearance of the scalar contraction of algebra generators, 
rVi = C^. This scalar product is known as the Casimir operator of the algebra. In the 
fundamental representation considered here, C^ ~ ^|l2- The Casimir clearly commutes 
with all group elements, thus (3.17) can be written as 

E.,{S)&{S)he[A] = -n^C^-f^ he,[A]heM] = -h^C^-fhelA]. (3.18) 

This expression will be useful below. 

On the other hand, if two consecutive fluxes act on one endpoint, say the target, we get 

E,iS)E,iS)he[A] = -ftV he[A]nT,. (3.19) 

From this result we immediately find that two flux operators do not commute, 

[E,{S),Ej{S)]K[A] = -h^^-" h,[A][T,,T,] = -hW^:3''he[A]Tk. (3.20) 

The actions (3.13) of the holonomy-flux algebra trivially extends to a generic basis ele- 
ment D^^\h). The action (3.13a) is unchanged, and in the right hand side of (3.13b) one 
simply has to replace t; by the generator Ji in the arbitrary irreducible j. Consequently, in 
(3.18) we have the Casimir C^ — —j{j + l)l2j+i on a generic irreducible representation, 

E,{SW{S)D'^^'>{h,) = n^j^jU + l)D'^^'>{he). (3.21) 

Finally, the action is extended by linearity over the whole Hkin- The remarkable fact is 
that this representation of the holonomy-flux algebra on Hkin is unique, as proved by Flei- 
schhack and Lewandowski, Okolow, Sahlmann, Thiemann [23]. This uniqueness result can 
be compared to the Von Neumann theorem in quantum mechanics on the uniqueness of the 
Schrodinger representation. It is well-known that the uniqueness does not extend to inter- 
acting field theories on flat spacetime. Remarkably, insisting on background-independence 
reintroduces such uniqueness also for a field theory. 

What we have accomplished with this construction is the definition of a well-behaved 
kinematical Hilbert space for general relativity. It carries a representation of the canonical 
Poisson algebra, and as a bonus, this representation is unique. Following Dirac, we now have 
a well-posed problem of reduction by the constraints : 

Hkin > T~l-kin ^ H-Diff > Hphys- (3.22) 
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3.1 Gauge-invariant Hilbert space 

The first step is to find the solutions of the quantum Gauss constraint. These are the states 
in Hkin that are SU{2) gauge invariant. These solutions define a new Hilbert space, that we 
call Tilin where we leave the subindex kin to keep in mind that there are still constraints to 
be solved before arriving to T-LpUys- The action of the Gauss constraint is easily represented 
in T-Lkin- In fact, recall that under gauge transformations 

he >K = Uohe = gs{e)heg^(ly (3.23) 

Similarly, in a generic irrep j we have 

D^^HK) -^ D^^HK) = D^'\gs(e)heg;(l)) = D^^Hg,f^,))D(^\K)D^^\g-l^). (3.24) 

From this it follows that gauge transformations act on the source and targets of the links, na- 
mely on the nodes of a graph. Imposing gauge-invariance then means requiring the cylindrical 
function to be invariant under action of the group at the nodes : 

fo{hi,...,hL) = fo{gs,higt,^^ , . . . , gs,^hLgtL~^) (3.25) 

This property can be easily implemented via group averaging : given an arbitrary / e Cylr, 
the function 

/o(/ii,...,/il) = / ]^dg„ /(5si^igti"\---,ffsi/iL5ti."^) (3-26) 

•' n 

clearly satisfies (3.25). 

Example : The theta graph. To give a constructive example, let us consider the following graph, 

.2 




A generic cylindrical function can be expressed in terms of the orthonormal basis using (3.12). Since 
the gauge transformations act only on the group elements, the gauge-invariant part is obtained looking 
at the gauge-invariant part of the product of Wigner matrices, 

Jinv\i^l^ • • • ^ 'f'S ) / ^ /mi ,7712 ,m3,ni ,n2,n3 -^TTiini I'^l J-^ 7712712 1*^2 J-l^ 7773713 1*^3 j 

L J inv 

Jl-"ii.n, 

Using the definition (3.26), the invariant part of the basis is 



^3) 

J inv 



= / 771i77127713QiQ20!3 ' /3l ^2 fe " 1 "2 »l3 "'^a i ^i ( "^ ) "'^02 fe ' '' "3 /Ss ' '^ •' ' 
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where 'Pmim2mr,aia2a3 IS the projector on the gauge invariant space, 

This projector can be written in terms of normalized Clebsch-Gordan coefficients, or Wigner's 3j-m 
symbols, as 



/' 



With this notation. 



s notation, 

L J inv 

h h J3\ fjl J2 J3\ fjl J2 J3\ fjl J2 J3\j^Ui) (, -.,^02) (, NpiOa) (, ^ 



ra\ 7712 7713/ \77i 712 713/-'--'- ' i^ 1. 

where j„ is a short-hand notation for the 3j-m symbols. Notice that these are the invariant tensors in the 
space of ® je of all the spins that enters in the node 77. Finally, we have 



, - VTTD'^^'Cft iTTi V f^i'^='« M' -^^ J3 \ Ml j2 T3 



ri m e 71 e 



^ ^ ;n,.2,.3 ■Q^Oe)(^^)-Q .^^ (3.28) 

with the new coefficients J^i J2 J3 including the sums over the magnetic numbers 777e, 7ie. 
The group averaging amounts to inserting on each node n the following projector, 

r^ fdgl[D^^^\g). (3.29) 

•^ e6Ti 

Here the integrand is an element in the tensor product of SU{2) irreducible representations, 

n^™l(Me(g)F(-'"=). (3.30) 

e e 

As such, it transforms non-trivially under gauge transformation and is in general reducible, 

(g)l^(^=)=0F(-'-). (3.31) 

e i 

Then, the integration in (3.29) selects the gauge invariant part of {^^ V^^" ^ , namely the singlet 
space V^'^\ if the latter exists. Since V is a projector, we can decompose it in terms of a basis 
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of V^^\ Denoting ia a vector ("ket") in this basis, a = 1, . . . , diuiV^^\ and ijj the dual 
("bra"), 

'P= Y. *"^«- (3.32) 

Q=l 

These invariants are called intertwine rs. For the case of a 3-valent node as in the above 
example, dini^'^^) = 1 and the unique intertwiner i is given by Wigner's 3j-m symbols 
(cf. (3.27)). More precisely in the case of a three-valent node the space 

yUi) YiJ2) ^ yijs) 

J inv 

is non-empty only when the following Clebsch-Gordan conditions hold, 

|J2-J3|<J1 <J2+J3. (3.33) 

For an n-valent node, the space V^^^ can have a larger dimension. To visualize the intertwi- 
ners, it is convenient to add first two irreps only, then the third, and so on. This gives rise to a 
decomposition over virtual links, which for n ~ A and n = 5 looks as follows : 





The virtual spins ki label the intertwiners. 



Brief proof that V is a projector. First we check that VP = V 

rr =fdg,dg2 n ^''='(51) n ^''^'(52) =fdgidg2 JJ D^'^\g,)D^'^\g2) = 

= fdgidg2 n D^'^\9i92) = f dg2 j dgi f] ^''^'(ffi) = V 
than that is left and right invariant 

l[D'^'^\gr)V= fdg2l[D^'^\g,)D^'^'>{g2)^ /d<72 H ^''^'(ffi.?^) = ^ 



(3.34) 
(3.35) 

(3.36) 



Vl[D 



Ue), 



dg2 n 0(^=)(52)Z?(^='(gi) = / dg2 n ^"^'(5251) = r (3.37) 



The facts that V acts only on the nodes of the graph that label the basis of Hkin and equa- 
tion (3.32) implies that the result of the action of V on elements of Hkin can be written as 
a linear combination of products of representation matrices D'^^\he) contracted with intert- 
winers, generalizing the result (3.28). The states labeled with a graph F, with an irreducible 
representation D^^^ [h) of spin-j of the holonomy h along each link, and with an element i of 
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the intertwiner space H,i = Inv[ (8) y^-'")] in each node, are called spin network states, and 

eGn 

are given by 

V'(r.j„,»„)[M = (E>D^'''>{he)(E}in. (3.38) 

e n 

Here the indices of the matrices and of the interwiners are hidden for simplicity of notation. 
Their contraction pattern can be easily reconstructed from the connectivity of the graph. 




To complete the discussion, let us show that imposing the gauge-invariance amounts to 
solving the Gauss law constraint, Giip = 0. Consider a gauge invariant node n, and a surface 
S centered in n of radius e. The action of the total flux operator through S* on n vanishes 
identically ; 

lim E{S) \n) = 0. (3.39) 

In fact, using (3.13b) at each link one notices that (3.39) produces the infinitesimal gauge 
transformation ga = I — ctiT^ G SU(2) at the node, and because the node is gauge invariant 
such action vanishes. 

Summarizing, spin network states (3.38) form a complete basis of the Hilbert space of 
solutions of the quantum Gauss law, 'H°^^. The structure of this space is nicely organized by 
the spin networks basis. As before, different graphs F select different orthogonal subspaces, 
thus 'H/Jjjj decomposes as a direct sum over spaces on a fixed graph. 



njO 



M. 



res 



(3.40) 



Furthermore, the Hilbert space on a fixed graph decomposes as a sum over intertwiner spaces. 



n'i = L2[SU{2f/SU{2)'^,d^iHaar] = ©j, (®nH„) 



(3.41) 



These two equations are the analogue in loop gravity of the Fock decomposition of the Hilbert 
space of a free field in Minkowski spacetime into a direct sum of n-particle states, and play 
an equally important fundamental role. 
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3.2 Geometric operators 

3.2.1 The area operator 

The simplest geometric operator that can be constructed in loop quantum gravity is the 
area operator. The area of a surface S can be given in terms of its normal Ua and the densitized 
triad E^, 

A{S) = / daidn2^ElE^^H^i, (3.42) 

The area at classical level. We start from the standard definition of area in terms of the metric, 



A{S)= j AarAa2\ldet[g^b^^] a,p = l,2 (3.43) 



usins the notation dix"' = ^ 

'^ OCT 

dx'^ dx 



det I dabTT-^jr-^ 1 = gab9cd \dix'^dix''d2x''d2x'^ - dix°'d2x''dix''d2x''^\ = (3.44) 



— gatgcd'^dix'^dix^ d2x'''d2X = 2g a[bg c]ddix'' dix d2x''d2X — (3.45) 
where we recognized that 



gg'^'ueUf (3.46) 



1 ef dx'' dx'' 

ga[bgc]d = -^Saceebdfgg rie = Eeab^-^^-^ (3.47) 

Using the definition of tetrad and densitized tetrad we have 

A{S) = I dcrido-2-y/e2efe/^nen/ = / A(TiAo2\J ElEi'^n^Uf (3.48) 

At the quantum level, we know from (3.13b) that the triad operators act as functional 
derivatives. The action of the scalar product of two triads operator was also studied, see 
equation (3.21), for the case of a surface intersected only once by the holonomy path. The 
case of a generic graph can be easily dealt with if we regularize the expression for the area in 
the following way. We introduce a decomposition of S" in A^ two-dimensional cells, and write 
the integral as the limit of a Riemann sum, 

A{S) ^ lim An{S), (3.49) 

where the Riemann sum can be expressed as 

N 

An{S)=J2VE,{Si)E^{Si). (3.50) 

7=1 

Here N is the number of cells, and Ei{Si) is the flux of Ei through the /-th cell. 
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Checking the limit. In the limit of infinitesimal cells we have that 

E,{Sl) ■: 



UaEid <y = EiUa ~ EfUaSl 

Sj J Si 



In that limit the definition of the area 

iV JV 



1=1 



(3.51) 



Am{,S) = Y, y/Er{Si)&{Si) ^ Y. ^ E-niSiE^^niSi = ^ 5/ ^ EfniE'^^ni = (3.52) 

1 7=1 1=1 

sjE^E^inant, (3.53) 



Accordingly, we define the area operator as 



AiS)^ lim An{S), 



(3.54) 



where in An{S) we simply replace the classical flux Ei{Si) by the operator Ei{Si). This 
operator now acts on a generic spin network state i/^r, where the graph F is generic and 
can intersect S many times. We already know that Ei{Si)E'^{Si) gives zero if Sj is not 
intersected by any link of the graph. Therefore once the decomposition is sufficiently fine 
so that each surface Sj is punctured once and only once, taking a further refinement has 
no consequences. Therefore, the limit amounts to simply sum the contributions of the finite 
number of punctures p of S* caused by the links of T. That is, 

N 



AiS) ^r = lim ^ JE,{SiWiSi) # = ^ hJj^jpijp + 1) Vr. (3.55) 



1=1 



pesur 



There are two key remarks to make to this formula : first of all, the spectrum of the area 
operator is completely known ^ and quantized : the area can only take up discrete values, 
with minimal excitation being proportional to the squared Planck length £p = hG, restoring 
Newton's constant. This result can be compared with other celebrated quantizations, such as 
the radii of electron's orbitals in atoms. 

Second, the operator has a diagonal action on spin networks. Therefore, spin network 
states are eigenstates of the area operator 

3.2.2 The volume operator 

Given a region i? C S classically we can define its volume as 



ViR) 






6:^x\ 



R 



Eabce'^^EfE^^El 



(3.56) 



4. In this expression, we assumed ttiat each puncture is caused by a link crossing the surface. However, it could 
also happen that S is puncture by a node of the graph. A closed expression for the ai'ea spectnjm is known also in 
this general case. See the literature for details [11, 13, 27]. 
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where the quantity in absolute value can be recognized as the determinant of the densitised 
triad, deti*;". 

Two distinct mathematically well-defined volume operators have been proposed in the 
literature. One is due to Rovelli and Smolin, and the other to Ashtekar and Lewandowski. We 
will refer to them as RS and AL respectively. Both of them act non-trivially only at the nodes 
of a spin network state. Let us begin reviewing the construction by RS. 

As we did for the area, we replace the integral over R by the limit of a Riemann sum. 
Specifically, we consider a partition of the region in cubic cells Cj so that R C U/C/, and the 
integral J„ d^x can be approximated from above by the sum J^i volume(C/). This partition 
allows us to rewrite (3.56) in terms of fluxes. In fact, consider the following integral, 

Wi = ^f d?ai f cfa^f d'<y3\s,jkEtiai)na{ai)E'^ia2)nbia2)E'k{a3W{a3)\ 
4o JaCi JdCi JdC'i 

In the continuum limit where we send the size of the cell e i-^ and we shrink the cell to a 
point X, we obtain 

Wi = -^e'''"=n„nbnedct£;f (x)e^ ~ dctEf{x)e^ ~ volumc2(C/). 
48 

Hence, we have 

V{R) = lim Y^ y/Wi. (3.57) 

For the sake of notation, let us subdivide each dCi into surfaces S*" such that dCi = 
\^aSJ ■ Then, we can write W7 as a sum of fluxes over three surfaces, and 



Finally, we can simply turn the classical fluxes to operators. 



(3.58) 




^(^)=]™oE./l^ E \e^,^E,{S'^)E,{S',)E^{S1) 



(3.59) 



This is the Rovelli-Smolin volume operator 

As for the area operator, one finds that there exists an "optimal" subdivision, after which 
the result stays unchanged with any further refinement, and so the limit can be safely taken. 
For the area, this consisted in the small surfaces being punctured only once at most. So- 
mething similar happens for the volume. The "optimal" partition is reached as follows. The 
nodes of F can fall only in the interior of cells, and a cell C/ contains at most one node. In 
case the cell contains no node, then we assume that it contains at most one link. Moreover, 
we assume that the partition of the surfaces dCi in cells SJ is refined so that links of F can 
intersect a cell S'j only in its interior and each cell S*" is punctured at most by one link. 

This said, let us now study the action of the operator The first thing to notice is that the 
presence of the epsilon tensor requires all three fluxes to be different : if two are the same, 
then their antisymmetric combination introduced by the epsilon vanishes. In particular, this 
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means that the volume does not act on links, since if no node is present, two of the Sf have to 
be the same. We thus obtain the important result that the volume operator acts only on nodes 
of the graph. 

Let us now focus on a single node, i.e. the /-th contribution to (3.59). We consider the 
cubic operator 



U^^y^ \e^,kE,{S'')E,{SP)Ek{S'^) 



48 ^ 

Let us restrict to gauge-invariant spin networks, where the Gauss law holds and each node is 
labeled by an intertwiner First of all, it is immediate to see that the action of (3.60) on a 3- 
valent node is zero. In fact, the Gauss law tells us that the sum of the fluxes through a surface 
around a gauge invariant node is zero. In the 3-valent case, only three S*" give non-zero fluxes, 
thus 

(e,{S'') + E,{S^) + E,{S^)) \i) = 0, (3.61) 

which implies 

E,{S-^) \i) = - (MS^) + E,iS^)) \i) . (3.62) 



Using this result in (3.60) we get zero because two identical fluxes always appear, 

e,jkMS'')E,{Sl')Ek{S^) K) = -e^jk (MS^) + E,{S^)) E,iS^)EkiS^) \i) = 0. 

Non-trivial contributions to the volume comes from nodes of valency 4 or higher Notice 
that these are the cases for which the intertwiner is not unique, but a genuine independent 
quantum number Thus the operator (3.60) probes precisely the degrees of freedom hidden 
in the intertwiners. Consider the 4-valent case. Of all the surface cells S", only four are 
punctured by the links. Hence, we have only four contributions to the Gauss law, which we 
can then use to eliminate one flux in favor of the remaining three, 

E,{S^) = -Ms') - Ms'-) - Ms"")- (3.63) 

The sum over the S" in (3.60) then reduces to the sole contributions from the four punctured 
surfaces. Using (3.63) these contributions are all equal. A simple combinatorial shows that 
there are 48 terms, all equal, thus 



U 



e^jkMs')EjiS^)EkiS'') ^h^\j^e,jkJlJf4\- (3-64) 



In the last step we have used the result (3 . 1 3b), which gives the action of the fluxes in terms of 
the generators J" in the spin ja representation. Notice that the orientation factors ± (between 
the edge and the surface) are irrelevant because of the modulus. 

This is a well-defined cubic operator, whose action spectrum is again discrete, with mini- 
mal excitation proportional to {hG)'^^^. The explicit formula is much more compUcated than 
that of the area. We refer the interested reader to the literature [29]. 

The alternative proposal by Ashtekar and Lewandowski has the extra property of being 
sensitive to the differential structure of the graph at the node. The subdivision of the region 
R into cubic cells is the same as above, the key difference is in how to choose the surfaces 
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to smear the triad filed. Instead of subdividing the boundary dCj into cells 5", we now 
consider for each C/ only three surfaces. We assign a local coordinate system .t°, and take 
the three surfaces to be any three surfaces S'^ inside the cube and orthogonal to each other. 
The Ashtekar and Lewandowski volume operator is defined as 



ni?) = limV 



y^,kea,.&{Sf)&{S])&{S'i; 



(3.65) 



Unlike (3.56), this expression is already written in terms of fluxes, so we do not need the 
manipulation in terms of Wi. We can directly use (3. 13b) to write 



1 



h^ 



^^e,jkeabcE\S^)E^ {S'})E''{Sj) - —e,,keabcKaKbK,r^JlJ^, (3.66) 

where again J* is the SU{2) generator in the spin ja representation, and ks is the relative 
orientation between the cell Sf and the edge puncturing it, i.e. 



1 if e lies above S 
K5(e) = < —1 if e lies below S" 
otherwise 



(3.67) 



The operator (3.66) strongly depends on the choice of the triple of surfaces Sf. In order 
to define a volume operator independent from the regulator, we average (3.66) on all the 
possible choices of orthogonal triples Sj (see the [27] for the details). The result is a sum 
over the nodes n of the graph. 



V{R) = ft3/2 J2 



riGr 



\ 



Kq 

48 



^ijk 



Y^ e(e,e',e") J'(e)JJ(e')J'=(e") 



(3.68) 



where e,e',e" runs over the set of edges passing through the node n and e is the orientation 
function, which equals if the tangent directions of the three edges are linearly dependent 
and ±1 if they are linearly independent and oriented positively or negative respect to the 
orientation of E. kq is an undetermined constant introduced by the averaging procedure. 

This alternative volume differs from (3.59) in two ways : firstly, it depends on an arbitrary 
constant, kq, unlike (3.59) which is unique. Secondly, the absolute value is outside the internal 
summation, unlike (3.59) where it is inside. This second differences is what makes the AL 
operator sensible to the differential structure of the graph ; in fact, the action of (3.68) vanishes 
on nodes whose edges lie on a plane, because in that case e = 0. 

On the other hand, both operators annihilate a 3-valent node, and coincide (up to a pro- 
portionality constant) in the 4-valent case, where they are both proportional to Sijk J} Jj J| . 

Summarizing, both volume operators act only on nodes of the graph. Their matrix ele- 
ments vanish between different intertwiner spaces, and since every intertwiner space is finite 
dimensional, their spectra are discrete with minimal excitations proportional to the Planck 
length cube £p. 
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Together with the discreteness of the area operator, these results show that in loop quan- 
tum gravity the space geometry is discrete at the Plack scale. Each spin network describes a 
quantum geometry, where each face dual to a Unk has an area proportional to the spin je, and 
each region around a node has a volume determined by the intertwiner i„ as well as the spins 
of the link sharing the node. 

It is important to stress that this is not a built-in discretization, as in lattice approaches to 
quantum gravity. It is a result of the quantum theory, similar to the quantization of the energy 
levels of an harmonic oscillator, or the radii of the atomic orbitals. 

Thanks to this fundamental discreteness, where the minimal geometric excitations are 
proportional to the Planck length, the theory is expected not to have ultraviolet divergences, 
and to resolve the problem of the classical singularities of general relativity at the big bang 
(e.g. [31]), or at the center of black holes (e.g. [32]). Addressing these issues of course re- 
quires taking into account the dynamics of the theory, to which we turn next. 



4 Loop quantum gravity : Dynamics 
4.1 Solutions of the diffeomorphisms constraint 

Spin network states '0(rje.j„)[^] ^^^ ^^ ^fchi' i-^- ^'V'(r.je,i„)[^] — 0- The next step in 
the Dirac program (3.22) is to implement the spatial diffeomorphisms, namely to find gauge- 
invariant states such that H°-il>[A\ = 0. 

To that end, consider a finite diffeomorphism (f). Its action on the holonomy, as in (2.46), 
naturally induces an operator cf) on the space of cylindrical functions. This operator maps 
Cylr to Cyl^oT, that is (jiipr = V'^or- Its action is well-defined and unitary, thanks to the 
fact that the Ashtekar-Lewandowski measure is diffeomorphism invariant. On the other hand, 
Cylr and CyZ^or are orthogonal Hilbert spaces, regardless of what the diffeomorphism is. 
This means that we can not define the action of an "infinitesimal" diffeomorphism diffeo- 
morphisms are all finite, from the perspective of cylindrical functions. Although this might 
appear as a restriction, it is not a real problem for the construction of T-Loiff- We can proceed 
by group averaging as we did for the Gauss constraint, and construct Hoiff from those states 
invariant under ^«/fe diffeomorphisms. 

However, there are a couple of subtle issues to take into account. The first one has to 
do with the existence of symmetries of the graphs. Namely, for each graph there are always 
some diffeomorphisms that act trivially on it, leaving it unchanged. Let us distinguish two 
cases : the diffeomorphisms that exchange the links among themselves without changing F, 
call them GSp following [30], and those that also preserve each link, and merely shuffle 
the points inside the link, call them TDiffp. The latter have to be taken out, because their 
infinite-dimensional trivial action would spoil the group averaging procedure. 

The second issue is that unlike imposing the Gauss law, imposing the invariance under dif- 
feomorphism will not result in a subspace of "H^j^j, since diffeomorphisms are a non-compact 
group. Think for instance of ip{x) e L2 [R, dx] required to be invariant under translations (a 
non-compact group) : the result is a constant function c, which is not in L2 [R, dx] since 
it is not integrable. It defines however a linear functional c : ip G L2 [R, dec] 1-^ C, since 
/ dxcipix) = c'i/'(0), the Fourier transform of ip evaluated in zero. Similarly, solutions to 
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the diffeomorphism constraint can be described in terms of linear functionab on 'HjJj^. Let 
us denote Ti.^*.,^ the space of linear functionals on "H^^n ■ Then, 77 € H°*„ is a diffeomorphic- 
invariant functional if 

r;(0V^) = r;(V.) Vt^gHL- (4.1) 

The space of such functionals is denoted Wj^^rr, and the Hilbert space Hoiff solution of 
the constraint is constructed by duality. The condition (4.1) should remind the reader of the 
analogue condition of gauge-invariance studied in the previous section. There, we were able 
to implement it via a simple group averaging procedure. In the same manner, we would like 
to define a projector Voiff on HDiff such that 



ii^mmff ^ {y^lVmffm 



E 

/.eDiff/TDiffr 



WIV^') 



(4.2) 



where the sum is over all the diffeomorphism mapping F into F' except those corresponding 
to the trivial ones TDiffp. See [30] for more details. 

The result of this procedure are spin network states defined on equivalence classes of 
graphs under dijfeomorphisms. These equivalence classes are called knots, see Figure 1 . The 
study of knots forms an elegant branch of mathematics. The diff-invariant Hilbert space of 
loop quantum gravity is spanned by knotted spin networks. 





Figure 1 - Left panel. A diffeomorphism can change the way a graph is embedded in E, but 
not the presence of knots within the graph. Right panel. The first few knots (without nodes), 
taken from Wikipedia. 



4.2 The Hamiltonian constraint 

Finally, we approach the last step of Dirac's program. On the space of knotted spin net- 
works Hoiff, we want to define the Hamiltonian constraint, and study its solutions. The 
classical scalar constraint is 



H{N) = 



d'xNe]^- 



EtE\ 



= //-^(iV) -2 (1 + 72) T(iV) 



F^,,-2{1+^^)KIkQ 



(4.3) 
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where we introduced the shorthand notation H-^{N) and T{N). As with the ADM Hamil- 
tonian constraint (1.19), this expression is non-linear, which anticipates difficulties to turn it 
into an operator. However, a trick due to Thiemann [33] allows us to rewrite (4.3) in a way 
amenable to quantization. Denoting V ^ J y^det{E) the volume of E, and 

K^JKEt, (4.4) 

we can use the classical brackets (2.34) to establish the following identities, 

K = ^{Al-K{E))=l.{Al,K}, (4.5) 



^ J fVk ^ ISA'' V\ (4 7) 

./Mm ""' ^t^-^^ (4./) 

Using these relations, 

H^'iN) = Jd'xNe-'-5,,F:,,{Ai,V}, (4.8) 

The advantage of this reformulation is that the non-linearity is mapped into Poisson bra- 
ckets. The next step is to rewrite these expressions in terms of holonomies and fluxes, so that 
we can turn them into operators. Notice that we already know the volume operator, and its 
spectrum can be explicitly computed. This is very promising towards the prospect of kno- 
wing the action of the Hamiltonian constraint. Next, the connection and curvature have to be 
written in terms of holonomies. This requires a regularization procedure. We describe it only 
for the first term, H^, and refer the reader to the literature [13] for T{N). 

The connection can be easily expressed in terms of holonomies. Writing explicitly (2.38), 
we have that for a path Cq of length e along the x" coordinate, h^^ [A] ~ 1 + eA^Ti + 0(e^ ), 
therefore we also have 

Kl {he^ ,V]^e {a;, V] + 0{^). (4.9) 

For the curvature, consider an infinitesimal triangular loop aab, lying on the plane ah, and 
with coordinate area e^. At lowest order in e, 

hc^. = 1 + l^'Kbr' + 0(e4), (4.10) 



thus 



2 



/ia„.-C -e'F>* + 0(e"). (4.11) 



At this point, we proceed as we did for the geometric operators in the previous section. We 
introduce a cellular decomposition of S, and regularize the integral as a Riemann sum over 
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the cells C/, 



H^ = lim y ef Se'^^^Tr {Fab {A^, V}) = 

/ 



(4.12) 
(4.13) 



This time is more convenient to specify the cellular decomposition in terms of a triangulation, 
namely a collection of tetrahedral cells. The loop aab can then be adapted to the triangular 
faces of this decomposition, as in the following figure. 




This expression can now be promoted to an operator in the quantum theory. 



H' 



I 



KL)h-' 



he^,V 



(4.14) 



This is a well-defined operator, whose action is explicitly known. It inherits the property of 
the volume operator of acting only on the nodes of the spin network. From the holonomies, 
it modifies the spin network by creating new links carrying spin 1/2 around the node, see the 
following figure. 




Finally, its amplitude depends on the details of the action of the volume operator See the 
literature [11, 13, 30] for details. 

Since the only dependence on e is on the position of the new link, in the Hilbert space 
T^Diff where the position of the link is urelevant, the limit can be safely taken without 
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affecting the result. Not only this operator is well-defined, also the Dirac algebra ( 1 .23) can be 
realized at the quantum level. To see this, we need first to note that the new "exceptional" links 
added by (4.14) carry zero volume and are thus invisible to a further action by (4.14). This 
crucial property requires the Ashtekar-Lewandowski version of the volume quantum operator, 
because the new nodes at the junction of the new link with the old ones are planar, and the 
AL volume operator vanishes on planar nodes. Thanks to this, H{Ni) and H{N2) commute 
on the space Ti-Diff of diffeomorphic-invariant spin networks. Hence for diffeomorphism 
invariant states ((/)| € T-Loiff, 

H(iVi),H(7V2)] IV') = = ((/.|H(iV[iViV2]) |V> , VV^ e Hfc,„. (4.15) 

This is the correct commutator algebra on-shell. Thanks to this non-trivial property, the quan- 
tization is said to be anomaly-free. 

Another key result is that an infinite number of states solutions of H are known : any 
graph without nodes is in the kernel of H^ and T. This result [35], together with the discre- 
teness of the spectra [34], started the whole interest in the loop quantum gravity approach. 
More in general, a formal solution of the Hamiltonian constraint will be a particular linear 
combination (generally infinite) of spin networks with an arbitrary number of exceptional 
links, whose coefficients depend on the details of the quantum operator and on the spins 
carried by the spin network. 



+ u; 



This construction defines a new spin network with a "dressed node" as the solution. However, 
this procedure is only formal, and no explicit solutions are known in general. 

In conclusion, we have a perfectly well-defined Hamiltonian constraint, whose action 
is explicitly known and finite. An infinite number of states solutions of H are known, and 
the Dirac algebra is anomaly-free on physical states. Compare this with the old-fashioned 
Wheeler-De Witt equation (1.37), which was badly ill-defined, and we see the full force of 
the use of the Ashtekar variables to quantize general relativity. 

Nevertheless, the program of quantization is far from being complete : the complete cha- 
racterization of Hphys nor the full spectrum of H are known. In spite of the successes of 
this approach, some limitations need to be stressed. In particular, although well-defined, the 
Hamiltonian constraint is plagued by a number of ambiguities : 

- We can change the spin of the exceptional edges to an arbitrary value j. This change 
have significant physical effects in loop quantum cosmology for example. 

- We can regularize the connection and the curvature with different paths and still make a 
consistent theory [30]. In particular, one can envisage a different construction in which 
the constraint acts on more than one node simultaneously, as advocated in [36]. 

- Alternative orderings can be explored. 
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4.3 Current approaches 

The effort to gain control over these ambiguities has flourished into two main lines of re- 
search. The first one is the idea of the Master Constraint. There one defines a unique constraint 
implementing simultaneously the diffeomorphisms and scalar constraints. The second one is 
the spin foam formalism. There one abandons the canonical approach, and seeks a func- 
tional integral description of transition amplitudes between spin network states. These two 
approaches have seen recent important developments, and are the current state of the art of 
the field. We will unfortunately not have time to review them here, but the interested reader 
can consider references [13, 37] for the first approach, and [38] for the second. 

To conclude, although the kinematics of loop quantum gravity is beautifully under control, 
the dynamics is still work in progress. We remark nonetheless that the key result on space dis- 
creteness, derived using kinematical states, is expected to hold also for physical states [11]. In 
the next and final Section, we want to go back to this discreteness and make some comments 
about it. 



5 On quantum geometry 

Let us review what we have learned in Section 3. The golden tools of loop quantum gra- 
vity are the spin network states. These form a basis in the kinematical Hilbert space, and 
diagonalize geometric operators. In particular, the quantum numbers carried by a spin net- 
work, {T,je,in), define a notion of quantum geometry,^ in the same way as the quantum 
number of an harmonic oscillator defines its quantum state. These quantum numbers can be 
compared with the kinematical metric gab, defining the classical geometry of space. We use 
the word kinematical to mean that gab is an arbitrary metric, not necessarily a solution of 
Einstein's equations, just like an arbitrary spin network spans the kinematical Hilbert space, 
not necessarily solving the diffeomorphisms and Hamiltonian constraints. The quantum geo- 
metry described by (F, je,in) is very different than a classical geometry gab- Specifically, it 
is largely insufficient to reconstruct a metric gab- We can highlight three key differences : 

(i) Quantized spectra : the spectra of geometric operators are discrete, as opposed to the 
continuum values of their classical counterparts. This is a standard situation in quantum 
mechanics, think for instance of the discretization of the energy levels of an harmonic 
oscillator. 

(ii) Non-commutativity : not all geometric operators commute among themselves. This is 
a consequence of the non-commutativity of the fluxes (3.20). This is also standard, like 
the incompatibility of position and momentum observables. 

(iii) Distributional nature : the states capture only a finite number of components of the 
original fields, that is their values along paths (for the connection) and surfaces (for the 
triad). This is reminiscent of what happens in lattice theories, where the continuum field 
theory is discretized on a fixed lattice and only a finite number of degrees of freedom 
are captured. 



5. Each face dual to a link has an area proportional to the spin je, and each region around a node has a volume 
determined by the intertwiner i„ as well as the spins of the link sharing the node. 
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In spite of these differences, if the theory is correct it must admit a semiclassical regime 
where a smooth geometry emerges. And furthermore, the dynamics of this smooth geometry 
should be given by general relativity, at least in some approximation. This is what we expect 
from any quantum theory : that the h ^^ limit is well-defined, and reproduces classical 
physics. The points (i-iii) show that this limit might not be trivial to obtain in loop quantum 
gravity. Although we started from genuine general relativity and quantum mechanics, without 
any exotic ingredient, we ended up with a description in terms of quantities, (F, je, in), which 
are very far from the original ones, gab- This is the problem of the semiclassical limit. 

In this final Section, we would like to comment on the first part of this problem, namely 
the emergence of a smooth geometry. This requires understanding the points (i-iii). The first 
point is easier to deal with : also the orbitals of the hydrogen atoms are quantized, placed at 
distances labeled by an integer n. The classical Keplerian behavior is recovered if we look at 
the large n limit. Similarly, continuum spectra are recovered in the large spin limit ji i—> oo. 
Points (ii) and (iii) are more subtle. The key to deal with them is the use of coherent states, 
namely linear superpositions of spin network states peaked on a smooth geometry. 

Recall that coherent states are peaked on a point in the classical phase space. For instance, 
the coherent state \z) for an harmonic oscillator is peaked on the initial position q — Re(2;) 
and momentum p = Im(2;). The phase space is the familiar cotangent bundle to the real line, 
{q,p) G T*R. Analogously, coherent states for loop quantum gravity are peaked on a point 
{A\{x) , Ef {x)) in the classical phase space, which defines an intrinsic (through the triad) 
and extrinsic (through the connection) 3-geometry. Such coherent states for loop quantum 
gravity were introduced in [39]. These states minimize the uncertainty of the flux operators, 
thus addressing (ii). In order to address (iii) and recover a smooth geometry everywhere on 
S, the coherent states have support over an infinite number of graphs. 

On the other hand, dealing with an infinite number of graphs is a formidable task, and 
for practical purposes, one often needs to rely on approximations. A convenient one is to 
allow only states living on a fixed finite graph F. The Hilbert space Hr provides a truncation 
of the theory, which may be sufficient to capture the physics of appropriate regimes. Within 
this truncation, one considers coherent states in T-Lr- However, in which sense one can assign 
a classical geometrical interpretation to these states ? They will be peaked on points in the 
phase space corresponding to "Hp, which consists of classical holonomies and fluxes on F. 
These quantities capture only ?l finite number of components of a continuum geometry. To be 
able to interpret these truncated coherent states in a physical sense, we need to use these data 
to approximate a continuum geometry. The problem is similar to approximating a continuous 
function f{x) if we are given a finite number of its values, /„ = f{xn)- Various interpolating 
procedures are common. We now describe a notion of interpolating geometry which emerges 
naturally. 

The associated gauge-invariant Hilbert space on a fixed graph is a sum over intertwiner 
spaces T-Ln = Inv[ ® V^^"''] associated to the nodes (see (3.41)), 

eGn 

^? = ©J, (®n^rO • (5.1) 

Just like L2[i?, dx] is the quantization of the classical phase space T*R, this Hilbert space 
is the quantization of a classical phase space ^, denote it Sy- The important result which is 



6. Up to singulai' points, but this is not important for the following. 
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relevant for us is that this phase space can be parametrized as follows [40], 

Sr= xT*S^ xSpin)- (5.2) 

I n 

Here T*S^ is the cotangent bundle to a circle, F the valency of the node n, and Sp is a 2(F — 
3)-dimensional phase space, introduced by Kapovich and Millson [41]. This parametrization 
mimics the decomposition (5.1). In particular, the quantization of T*S^ gives the quantum 
number ji associated to each link, and the quantization of 5i?(„) gives the F — 3 intertwiners 
associated to an F-valent node. 

The parametrization defines a notion of interpolating geometry associated to a cellular 
decomposition of S dual to F, called twisted geometry [40]. The key to the geometric inter- 
pretation of (5.2) is the node space 5f-(„) . Fix the spins ji on the links connected to the node, 
and consider F unit vectors ni in M'^ constrained to satisfy the following closure condition, 

F 

C = Y,Jini = 0- (5-3) 

1=1 

Because of the closure, the vectors jini span a polygon in M'^, see Figure 2. The space of all 



J2"'2 




Figure 2 - A polygon in R"^, with F sides of fixed lengths. Because of the closure condition, 
not all the components of the normals are independent, but only 2F — 3. If we further consider 
the polygon up to rotations, that is the space Sf, this is only specified by 2F— 3— 3 = 2(F— 3) 
numbers. These can be taken to be the lengths of F— 3 diagonals and their conjugated dihedral 
angles (which measure the way the polygon hinges on each diagonal). 

polygons with fixed ji up to SO(3) rotations forms the phase space Sp, called by Kapovich 
and Millson the space of shapes of the polygon. 

This means that we can interpret the intertwiners as "quantized polygons". Such inter- 
pretation is appealing, but not particularly useful for our purposes. More promising would be 
a description in terms of polyhedra, rather than polygons : Polyhedra can be glued together 
to approximate a smooth manifold, and their geometry will induce a discrete metric of some 
sort. 

The good news [43] is that this is precisely what happens ! In fact, each (non-coplanar) 
configuration of vectors ni describes also a unique polyhedron in M."^, with ji as the areas of 
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Figure 3 - A generic polyhedron with 100 faces. Notice that most faces are hexagons. 



the faces, and n; as their unit normals. More precisely, a convex bounded polyhedron. The 
explicit reconstruction of the polyhedron geometry from holonomies and fluxes is studied in 
[43]. Therefore, up to the degenerate configurations corresponding to coplanar normals, we 
can visualize Sp as a polyhedron with F faces. 

Going back to (5.2), the T*S^ space on the links carries information needed to define the 
parallel transport between adjacent polyhedra. This unique correspondence between Sp in 
(5.2) and polyhedra means that each (non-degenerate) classical holonomy-flux configuration 
on a fixed graph can be visualized as a collection of adjacent polyhedra, each one with its 
own frame, and with a notion of parallel transport between any two polyhedra. Since twisted 
geometries parametrize the classical phase space on a fixed graph, coherent states on a fixed 
graph are peaked on the approximated smooth geometry described by a twisted geometry, i.e. 
a collection of polyhedra. 

Twisted geometries have a peculiar characteristic which justifies their name : they define 
a metric which is locally flat, but discontinuous. To understand this point, consider the link 
shared by two nodes. Its dual face has area proportional to ji . However, its shape is determined 
independently by the data around each node (the normals and the other areas), thus generic 
configurations will give different shapes. In other words, the reconstruction of two polyhedra 
from holonomies and fluxes does not guarantee that the shapes of shared faces match. Hence, 
the metric of twisted geometries is discontinuous in the sense that the shape of a face "jumps" 
when going from one polyhedra to the next. See left panel of Figure 4. 

Notice that one can also consider a special set of configurations for which the shapes 
match, see right panel of Figure 4. This is a subset of St which corresponds to piecewise 
flat and continuous metrics. For the special case in which all the polyhedra are tetrahedra, 
this is the set-up of Regge calculus, and those holonomies and fluxes describe a 3d Regge 
geometry. The matching conditions in this case were studied in [44]. For an arbitrary graph, 
the matching conditions are studied in [43], and the result would be Regge calculus on generic 
cellular decompositions. 

The relation between twisted geometry and Regge calculus implies that holonomies and 
fluxes carry more information than the phase space of Regge calculus. We stress that this 
is not in contradiction with the fact that the Regge variables and the LQG variables on a 
fixed graph both provide a truncation of general relativity ; simply, they define two distinct 
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Figure 4 - Two high-valency nodes and the link connecting them. Left panel : In general, 
the two adjacent polyhedra, defined by holonomies and fluxes on the nodes, don't glue well. 
Even if the area of a polygonal face is the same because there is a unique spin associated 
with the link, the shapes will be different in general. Right panel : in order for the shapes to 
match, one needs to impose appropriate conditions on the polyhedra such as the matching of 
the normals to the edge in the plane of the face. These conditions, studied in [43], affect the 
global shape of the polyhedron. 



truncations of the full theory. See [45] for a discussion of these aspects. 

Summarizing, when we truncate the theory to a single graph, we capture only a finite num- 
ber of degrees of freedom of the geometry. This finite amount of information can be suitably 
interpolated to give an approximate description of a smooth geometry. The reparametriza- 
tion in terms of twisted geometries achieves a prescription for this interpolation procedure : 
Starting from holonomies and fluxes on a graph, we can assign to them a specific twisted 
geometry, that is, a bunch of polyhedra stuck together 

These results offers on the one hand a new way to visualize the quantum geometry of 
loop quantum gravity in terms of fuzzy polyhedra. On the other hand, they open the door to 
new beautiful mathematical ingredients, from the geometry of polyhedra [42, 43], to twistors 
[46], to new emerging symmetries [47]. 

These intriguing aspects, together with the striking recent developments in the study of 
the dynamics and in the applications to cosmology and the physics of black holes, make of 
loop quantum gravity an excellent arena for a student interested in the problems of physics at 
the Planck scale. 
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